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Casimir-Lifshitz interaction emerging from relative movement of layers in stratified dielectric 
media (e.g., non- uniformly moving fluids) is considered. It is shown that such movement may result 
in a repulsive Casimir-Lifshitz force exerted on the layers, with the simplest possible structure 
consisting of three adjacent layers of the same dielectric medium, where the middle one is stationary 
and the other two are sliding along a direction parallel to the interfaces of the layers. 
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I. INTRODUCTION 

In this paper we consider Casimir-Lifshitz forces [1, 2] in layered moving media. Our interest to this problem was 
initiated by a recent discussion on the friction forces that may or may not appear due to quantum-electromagnetic 
fluctuations in systems involving moving dielectric slabs [3-8]. In this paper, however, we will concentrate on another 
interesting theoretical issue which, to the best of our knowledge, has not been addressed so far: on the possibiUty 
of having repulsive Casimir-Lifshitz forces in moving dielectrics. The so-called Casimir repulsion is known to appear 
between electrically and magnetically polarizable objects in vacuum [9-12], or between dielectric objects of different 
permittivity that arc immersed in a dielectric fluid of an intermediate permittivity [13- 15]. Very recently, ultralong- 
range repulsive forces in piston configurations involving cut metallic nanorods have been reported [16]. There have 
been also attempts on achieving repulsion or "quantum levitation" with the use of other metamaterials [17-24]. 
However, recently it has been shown that the force between mctal-diclcctric mctamatcrial slabs in vacuum is always 
attractive [25-27]. The symmetry considerations also impose restrictions on the sign of the Casimir force [28, 29]. 

In this paper we are going to consider the case in which the force appears only as the result of relative movement of 
dielectric layers. In contrast to the Casimir friction studies, we are interested in the force component perpendicular 
to the direction of the movement. The main idea of this work is to consider a system which is initially balanced, i.e., 
when there is no movement there are no fluctuation-induced forces. One example of a system with such property is a 
uniform medium, say a fluid, which is initially at rest. There is, however, a possibility that when separate layers of a 
fluid begin to slide one with respect to another the balance is destroyed and there appears a noncompensated attractive 
or repulsive interaction between the sliding layers. It should be well understood at this point that the situation that 
we consider in this paper differs principally from the previously studied case of moving dielectric slabs separated by a 
vacuum [3, 4]. In the latter case, with an appropriate Lorentz transformation for the electromagnetic field, one may 
always reduce a problem involving a moving slab of an isotropic dielectric in vacuum, to an equivalent problem with a 
stationary slab of the same isotropic dielectric in vacuum. This is possible because under a Lorentz boost the vacuum 
"background" remains itself. Quite differently, in this paper we study the Casimir-Lifshitz interactions that appear 
in non-uniformly moving matter. Applying Lorentz transformations in this case results in a more difficult problem 
involving layers of anisotropic and nonreciprocal media. 

Therefore, we are going to approach this problem without resorting to an assumption that the available theories 
[1, 2, 30, 31] of the Casimir-Lifshitz forces in dielectrics are also applicable in the case of moving media. Instead, 
we quantize the electromagnetic field in moving matter and derive a relation for the zero-point energy from the 
first principles. This is required because moving media arc not invariant under time reversal and the traditional 
quantization scheme based on a modal expansion in a large box is not applicable (at least, without significant 
modifications). In fact, in this work we develop an alternative quantization approach that allows to reuse many 
of the results of the classic treatment of such nonreciprocal media. Nevertheless, the results of our method fully agree 
with the phenomenological quantization schemes developed by other authors [32-34]. 

The nonreciprocity considered in this paper is twofold: it may either be a result of relativistic movements of material 
fiuids or it may manifest itself in uniaxial bianisotropic metamaterials the constitutive relations of which include a 
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term that is responsible for nonreciprocal magnetoelectric coupling. Such mctamaterials have been theoretically 
known for a long time [35-37]; certain practical realizations have been proposed as well [38]. Some authors do not 
make a clear distinction between the real moving media and their metamaterial counterparts, taking for granted 
that the two types can be described with the constitutive relations of the same form. This is, however, not entirely 
true. Although applying the Lorentz transformations to the Maxwell equations written for a moving dielectric results 
(in the laboratory frame) in bianisotropic material relations with nonreciprocal magnetoelectric coupling, such a 
transformation may not always lead to spatially local constitutive relations. Indeed, the Lorentz transformation 
intermixes the spatial coordinates with time, therefore, a medium which is nonlocal in time in one of the reference 
frames (i.e., a dispersive dielectric in its proper frame) becomes nonlocal in both space and time in another reference 
frame. Thus, a moving dispersive dielectric may be described (in the laboratory frame) with the equivalent spatially 
local bianisotropic material relations only in a limited frequency range where the dispersion is negligible. 

Therefore, in this work the emphasis is mostly on weakly dispersive moving magnetodielectrics for which one may 
assume that s{uj) and /i(w) are practically constant and real in a wide range of frequencies. This simplification, 
however, is not crucial for the main theoretical prediction of this paper, namely, the existence of repulsive Casimir- 
Lifshitz forces in layered moving media. This can be seen from the known fact (see, e.g. [20]) that the range of 
frequencies that make the dominant contribution to the Casimir energy in a pair of material layers separated by the 
distance d is limited by [wmaxj ~ '^irVp^/d, where z;ph is the phase velocity in the background material. Thus, if Wmax 
is set to the upper boundary of the region of low dispersion of a medium, then the theory developed in this paper 
will apply at separations d > 27rwph/|wmax|- As there exist real materials with low dispersion and loss up to, at least, 
the ultraviolet band, the applicability range of our theory may start at hundreds of nanometers. A straightforward 
generalization of the theory to the dispersive case is outlined in one of the appendices. 

The paper is organized as follows. In Section II we solve classically for the eigenwaves in a moving nondispersive 
medium and discuss their properties. In Section III we derive an expression for the Hamiltonian of the free elec- 
tromagnetic field in a moving medium and prove an orthogonality relation that holds for the eigenmodes in such 
a nonreciprocal medium. In Section IV we quantize the (macroscopic) electromagnetic fields in a moving medium 
and express the Hamiltonian of the electromagnetic field in terms of the creation and annihilation operators of a 
bosonic field. In Section V we obtain an expression for the zero-point energy and its regular part that represents 
the Casimir-Lifshitz interaction energy. In Section VI we solve for the Casimir-Lifshitz force in layered moving me- 
dia. In Section VII we present and discuss some numerical results that clearly demonstrate existence of repulsive 
Casimir-Lifshitz forces in such media. 



II. ELECTROMAGNETIC WAVES IN A MOVING MEDIUM 



We consider a uniaxial medium (the axis is along Zq) which is characterized by material relations of the following 
form [in this section we work in the frequency domain; the time dependence is of the form exp(— iwf)]: 

D = f • E-hazo X H, (1) 
B = p • H - azo X E, (2) 

where e = e^It + ezozo and /tZ = fi^It + H'z.q'Lq are the dyadic permittivity and the permeability, respectively, with 
being the unity dyadic in the plane transversal to zq, and a is the parameter of magnetoelectric coupling. Notice that 
due to the choice of signs in (l)-(2) this coupling in nonreciprocal. 

Such a medium can be envisioned either as a metamaterial with nonreciprocal bianisotropic inclusions, or as a 
an effective medium resulting from application of the Lorentz transformations to the electromagnetic fields in a 
magnetodiclectric moving with certain velocity v along the 2;-axis. In the latter case, the material parameters as seen 
in the stationary frame satisfy (see, e.g., [39]) 



'1 - n2/3 
/3 n2 - 1 
c 1 — n' 



/^t = M^ 1;^, (4) 

2«2' (^) 



where e and ji are the permittivity and the permeability in the comoving frame, c = 1 / y^eoMo is the speed of light in 
vacuum, /3 = v/c, and = £/x/(eo/xo). The material parameters are assumed nondispersive and lossless in (3)-(5), 
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but, in fact, these relations may be also generalized for dispersive moving media if plane waves are considered (this is 
further discussed in Appendix C). 

It should be noted that when these transformations are applied to a medium with = 1. they result in a = 
and the old values of the permittivity and permeability, independently of the velocity v. Thus, due to (3)-(5), a 
vacuum appears as a "medium" with properties invariant with respect to relative motion, while media with nontrivial 
refractive index are seen differently in different inertial frames of reference. 

The Maxwell equations for the fields in a moving medium can be written as 

iojji • H = Vt X E + {ioja + dz)zo x E, (6) 
-iuE • E = Vt X H + {ioja + dz)zo x H, (7) 

where Vt = A • V and dz = d/dz. Seeking for plane wave solutions of (6)-(7), it is possible to reduce Eqs. (6)-(7) to 

uj^etiit - {k. + wa)2 - ^fc^l ^ ^^^^^^ 



uj StfH - {kz + uja) k^ 



Hz=0, (TE,), (9) 



where k = kt + fc^ZoZo, kt = /t • k, is the wave vector of a plane wave, and the two equations (8) and (9) are 
for two independent polarizations: the transverse magnetic polarization with respect to the 2;-axis (TM^), for which 
Hz = 0, and the transverse electric polarization (TE^), for which Ez = 0. The transversal components of the electric 
and magnetic fields (with respect to the 2;-axis) in both TM^, and TE^. polarizations can be expressed through the 
^-components of the fields: 



xjTM. _ ^£(kt X zq) pTM. _ <kz + cja)kt ^ 

ilt - p ^z, Et Ez, (J-M^j, (lUj 



TE, a;M(kt x zp) te, Kkz + ^a)kt „ . , . 

* = ^2 = ^z, (TE^)- (11) 



The electric displacement D and the magnetic induction B in the same modes can be found with the help of the 
material relations (l)-(2) and the relations (lO)-(ll): 



TM, _ £ ((£tMt - - akz){kt x Zq) 
sk^'kt 



^' * (TM,), (12) 



^TE, _ M((£tMt - a^ )u; - afc^)(kt x Zq) ^ ^te. 



ak k ^^^^ 

BTE^ _ M^z^t jj dTE^ _ 

t — I3~^z, -t>z — M-f^zj 

"^t 

An interesting property of the D and B vectors in a moving medium is that despite the fact that the medium is 
anisotropic the three vectors k, D, and B are mutually orthogonal in each of the TM^ and TE^ modes. 

In the nondispersive case the dispersion equations (8) (9) arc quadratic with respect to the frequency and can be 
easily solved. As follows from (3)-(5) and (8)-(9) the equations are the same for both TM^ and TE^ modes. The 
roots of the dispersion equations are given by 

c^i,2(k) ^ Pjn^ - l)kz ± V(l - P^n^ - P^)kl + 71^(1 - /32)fc2] 

c n2-/32 • ^^^> 

The expression under the square root is nonnegative because >1 and < 1. In the case when P'^v? < 1 (i.e., 
when the velocity is below the threshold of the Chcrenkov effect) only a single solution of the dispersion equation is 
nonnegative, namely, the one with the plus sign in (14). 

When 0^11? > 1 there may exist zero, one, or two nonnegative roots of the dispersion equation, depending on the 
wave vector. Without any loss of generality we may assume v > and, thus. 6 > 0. Then, the roots of the dispersion 

equation are both negative (positive) if fc^ < (fcz > 0) and \kz/kt\ > \J fpn^l^i - When |fc^/fct| < y^ ^2~f_^f , there 
are two roots of opposite signs. The boundary between these regions defines the Cherenkov cone as seen from the 
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stationary frame. In the comoving frame (i.e., in the frame in which the medium is at rest), the same cone is seen as 
having the half-angle 6 such that tan^^ = . =, which is a well-known result. 

The nonrcciprocity of the material relations (1) (2) results in an obvious property of Eqs. (8)-(9): these equations 
are not invariant with respect to the change of sign of uj. However, the time-harmonic fields that are the solutions 
of (6)-(7) must satisfy the reality condition F_t^(x) = F* (x), where F represents either E or H and the symbol * 
denotes complex conjugation. Thus, their spatial Fourier transforms, F^{k) = / (fixFi^{y:)e~^^'^ , that represent the 
complex amplitudes of the respective plane waves, are such that F_-j(— k) = F^(k). It is immediately seen that the 
equations (8)-(13) are invariant under such a transformation that changes the signs of u and k simultaneously. In 
addition to this, the dispersion equations are also invariant with respect to a simultaneous change of signs of co and 
kz, as follows from (8)-(9). 

The instantaneous fields in a given polarization, F(x, t), can be written as a superposition of the plane wave solutions 
of (8) or (9): 



F(x,t) = J ^^F„^(k)e^C— *), 



(2^)^ V 



(15) 



where the index p = 1,2 labels the roots ujp — a;p(k) [Eq. (14)] for a given k, and Fi^^{]i) represent the complex 
amplitudes of the waves that belong to the two different branches of (14). 

The reality condition F_tj(— k) = F* (k) allows to rewrite (15) as follows. We notice that the two branches of (14) 
are such that wi(k) = — W2(— k), and W2(k) = — wi(— k). Hence, by replacing k with — k in one of the addends of the 
sum in (15), Eq. (15) can be written in the following equivalent form where only a single branch occurs explicitly: 



F(x,i) 



(2^ 



F<,(k)e*('^-^-"*)+F:(k)e 



—i{k.-x—Ljt) 



(16) 



Any branch may be chosen; for the following we select the branch with the plus sign in front of the square root in 
(14). 



III. THE HAMILTONIAN OF THE FREE ELECTROMAGNETIC FIELD 



Classically, the Hamiltonian of the free electromagnetic field in a moving medium can be obtained by considering 
the Maxwell equations written for instantaneous fields: 



dtB = -V X E, 
atD = V X H, 

where dt = d/dt. Performing the standard steps on derivation of the Poynting theorem, we write 
V • (E X H) = (V X E) • H - (V X H) • E = -{dtB) ■ H - (9*0) • E. 



(17) 
(18) 



(19) 



Next, we use the material relations (l)-(2) to express D and B in terms of E and H and, after recollecting the terms 
on the right-hand side with some trivial vector algebra, we obtain 



V • (E X H) = -dt 



B H D E 



(20) 
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i.e., the same final result as in a stationary medium. Thus, the Hamiltonian is (the same expression was used in [33]) 



B H D E 

2 ^ 2 



d^k 



B(k, t) ■ H(-k, t) D(-k, t) ■ E(k, t) 
2 ^ 2 



(21) 



where F(k, t) (with F representing any of the fields) are the time-dependent spatial Fourier transforms defined by 
(15): 



F(k,f)=^F,^(k)e- 



(22) 



Now we substitute the above representation into (21) and obtain 
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1 C fl^\c 



1 f d^k 

2] Wf'^ 



+ 



[B^^(k).H:^(k) + D:^(k).E„^(k) 

2 j (^E[B-.(k)-H:^(k)+D:^(k).E„^(k)]e-^(---=)*. (23) 



Since the dispersion equations are the same for both TE^ and TM^ modes, the field vectors that appear in (23) may 
be regarded as arbitrary hnear combinations of the fields of these two main polarizations. 

In an isolated conservative system the expression (23) represents the total electromagnetic energy that remains 
constant when the system evolves with time. Therefore, the last integral term in (23) that explicitly depends on time 
must vanish. Using the Maxwell equations written for plane waves, this term can be expressed as 



2 



^ E [B^^ (k) . (k) + (k) . E^^ (k)] e- 



1 /■ d^k 

27 Wrh 



^k. [D:^(k) X B^^(k)]e-*(-''-=)*, (24) 



where fc = |k|. It is immediately seen that in a reciprocal medium this integral vanishes for arbitrary Fourier 

transformed fields, because in such a medium LOs(k) = LUs{—k) = — a;p(k). 

In the moving medium, however, the situation is more complicated. Consider, for example, the case when at t = 
the electromagnetic field forms a pulse composed of waves with the wave vectors concentrated around k = +ko and 
k = — ko. This situation corresponds to defining an initial condition for the fields in terms of an oscillating function 
(oscillating in space!) with a smoothly varying amplitude vanishing at infinity. Then, in this pulse there are waves 
with frequencies concentrated around co = ±cji(ko) and uj = ±W2(ko), whereas wi(ko) — W2(ko). Let us look closer 
at the term (24) in this case. We may get rid of the integration around ±ko in (24) because the spectral width of the 
pulse is assumed to be small. Dropping an insignificant constant factor we obtain 



2 ^ ^ 

^i". k=±ko 



p,s ■ 

p^8 



^k.[D:^(k)xB„^(k)]e- 

2 ^ "^^^^^ • ^-^(ko) X B„^(ko) +D„^(ko) X B:,(ko)]e-'("--»)* = 

= ^^^^^Re{ko • [D:^(ko) X B,,(ko) +D,,(ko) x B:^(ko)]e-^("^-"^)*} . (25) 

The only possibility to make this term independent of time is to have its amplitude vanishing: 

ko • [D:^ (ko) X B^, (ko) + D^, (ko) X B:^ (ko)] = 0. (26) 

It can be verified by direct substitution that this condition holds for both TE^ and TM^; modes (and also for any 
linear combination of them). Because the same transformation that we have done above could be applied directly to 
the integrand of (24), we have proven that the term (24) vanishes in general. Physically, Eq. (26) has the meaning of 
an orthogonality condition for the modes with the wave vector ko and the frequencies a;i,2(ko) in a moving medium. 
Thus, we have proven that the Hamiltonian in a lossless non-dispersive moving medium can be written as 



^ = ^ / E • H:,(k) + D:^(k) . E„^(k); = 

= / -0^sY.'i^-^^l^pO-)^^^.(^)\= j (0gk.[D:(k)xBUk)]+c.c., (27) 
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where in the last equaUty only a single branch occurs as in (16), and "c.c." denotes the complex conjugate of the first 
term. The only difference of (27) from the same expression for a reciprocal magnetodielectric is in that a;(— k) 7^ w(k). 

We can split the total electric and magnetic fields in the last expression into the components corresponding to the 
TM^ and TE^ polarizations. To do this we notice that the electric displacement vector of the TM^ (TE^) mode and 
the magnetic induction vector of the TE^ (TM^) mode are collinear. Thus, the cross terms in the vector product 
D* (k) X B(^(k) do not contribute to the Hamiltonian (27). Therefore, we can write 

^ = / (C^p(k-[D:(k)xB„(k)]™.+k.[D*(k)xB„(k)]™^)+c.c., (28) 

where the brackets [. . .] [. . .]'^^^ denote the separate contributions of the respective modes. 

The relations (27) and (28) have a clear physical meaning. Indeed, the term (k/fc) • [D* (k) x Bj^. (k)] corresponds 
to the momentum of a plane wave in the moving medium. From the other hand, the energy w and the momentum p 
of a plane wave arc related hy w = {uj/k)p. Therefore, Eqs. (27)-(28) may be understood as a summation over the 
energies of all possible plane waves. 



IV. QUANTIZATION OF THE ELECTROMAGNETIC FIELD IN A MOVING MEDIUM 

The quantization of electromagnetic field in a moving medium is a well-established subject (at least, in the nondis- 
persive case) and can be performed within different frameworks: (i) with the covariant Lagrangian formalism of 
Ref. [32] , (ii) with the Heisenberg formalism of Ref . [33] , (iii) with the Green tensor-based formalism of Ref. [34] . All 
these approaches agree and lead in effect to the so-called canonical quantization of the electromagnetic field. 

Perhaps, the most intuitive approach is the one based on Heisenberg formalism. Under this approach one starts 
with a classical expression for the Hamiltonian H in terms of the instantaneous fields as in (21). The field variables 
as functions of the position and time that appear in the Hamiltonian are promoted to Hermitian operators that 
satisfy certain commutation relations. The commutation relations must be such that the equations of motion in the 
Heisenberg formalism (here and in what follows the square brackets [•,•] denote the commutator of two operators: 
[A, B]=AB- BA) 



(29) 
(30) 



result in a system of partial differential equations identical in form with the classic Maxwell equations. Such an 
equivalence exists because the classic electromagnetic theory can be thought of as a theory of quantum states of light 
with a very large number of photons (which are bosons) in each state. Although not required in the time-harmonic 

regime, for an arbitrary time evolution the above (sourccless) equations must be complemented by V • D = V • B = 0. 

In [33] it was shown that the required equal-time commutation relations can be written in terms the Cartesian 
components of B and D as 



[A(x,<),Bj (x',i)] = -iheijkdk5{yi- -x.'), 



(31) 



where Sijk is the Levi-Civita tensor, dk = d/dxk, and 6{x) is the three-dimensional Dirac delta function. Here and in 
what follows we use Einstein's notation in which a summation over repeating indices is assumed. It is also assumed 
that all components of D commute in between themselves, as do the components of B. 

From the commutation relation (31) it is seen that the noncommuting components of the field operators B and D are 
mutually orthogonal. Let us show that indeed such commutation relations lead to the Maxwell equations (17)-(18). 
First, we rewrite the Hamiltonian of the electromagnetic field (21) in terms of only B and D: 



■H= d'yi 



B 



B D £ D 

— + 1 



XZo • (B X D - D X B) 



(32) 



where ?? = (et/xt - a^) ^eth + IJ- ^zqZo, £. = (etMt - a^) VtA 
of the material relations (l)-(2) transformed to the form 



-e ^zqZo, and x = a(etAtt 



are the parameters 



E = £ • D - XZo X B, 
H = ^ • B -I- XZo X D. 



(33) 
(34) 
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The symmetry of ^ and rj and the form of the last addend under the integral (32) provide that the Hamiltonian is a 
self-adjoint (Hermitian) operator: T-L^ =T-L (here and in what follows the symbol ^ denotes Hermitian conjugation). 
Then, calculating, for example, the commutator of B and H we find 

[Bi,U] = ^J rf^x[Bi(x'), (^«/3^a(x)I?/3(x) - X£za0{Ba{^)D^{K) + I?^(x)S„(x))] = 

d3x(^„^([Bi(x'),£>„(x)]£»;3 + D^[Bi{^'),D0{x)]) - 2x£.a/3S„(x)[Bi(x'),£';3(x)]) = 

ih f 

= y / d-^xdkd{x - x') [ia/}{£aikD/}{x) + £^ifcD„(x)) - 2x£za/3e/3iftS„(x)] = 

= -ih j d^x6{x - x')Sikadk [^a/3-D^(x) - X£azpBp{x)] = -ihSikadkEa, (35) 

which is the same as [B,?^] = —iHV x E. In the derivation we used the fact that = ^pa- In a similar manner one 
obtains [T>,n] = ih"^ x H. 

The standard way to proceed after this step is to make a transition into the momentum space by expressing D(x) 
and B(x) in terms of a pair of conjugate canonical variables P(k) and Q(k). One then diagonalizes the Hamiltonian 

written in terms of P(k) and Q(k) by introducing the creation and annihilation operators. We would like, however, 
to move along another way that will allow us to reuse many of the results of the classic theory considered in the 
previous sections. 

To make a connection with the frequency domain treatment of Section II we look for the solutions of the Maxwell 
equations (written for the quantum vector field operators!) that have the form (here F represents any field vector) 

/rP\< 
^F^(k)e^[''— -C')*!, (36) 

where the operators Fa,(k) can be understood as the (time and position-independent) wave amplitude operators. 
The reality condition requires F(x, t) to be an Hermitian operator, thus F_(^(— k) = Fj,(k). When such a form 
is substituted into the Maxwell equations, one can reduce these equations to (8)-(13) with all the field variables 
promoted to wave amplitude operators. 

As the wave amplitude operators are assumed non-trivial, the frequency w(k) in (36) is found by solving a dispersion 
equation which is identical to the classic one. Thus, there are two dispersion branches a;p(k) = —ujs{—i^), s ^ p, and, 
analogously to (15) and (16), we can write 

F(x,t) = / [F„(k)e'('^— +Ft(k)e-'(''— , (37) 

where only a single branch w(k) appears explicitly [as before, we select the branch with the positive square root in 
(14)]. 

With enough care, the results of Section III may be also promoted to operators, provided that they are written in a 
form that satisfies the reality condition for the Hamiltonian: T-L^ = T-L. Thus, in the operator form the Hamiltonian (27) 
becomes 



1 r rjilr^ 

^=27 (^pk.[Dt(k)xB„(k)-B^(k)xDi(k)]+h.c.= 

= 2 y (^^[[kxDt,(k)].B„(k)+B„(k).[kxDt,(k)]]+h.c., (38) 



where "h.c." stands for the Hermitian conjugate of the first term. 

The representation of the Hamiltonian in terms of D and B is useful because in both TM^ and TE^ modes in a 
moving medium the three vectors k, D, and B are mutually orthogonal, as has been found in Section II. Therefore, 
in each mode separately the vectors B and k x D are collinear, while the same vectors corresponding to the two 
different modes are mutually orthogonal. Thus, with the help of Eqs. (12)-(13) we may express the vectors B„(k) 
and k x Dt^(k) as 



B„« = \/5^(.o.«e..^e,|, (3.) 
k X D^^(k) = kJ (^^^^ ®i 7oCv'etMta2(k) 62 ) , (40) 
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where Jq = [(et/^t — a^)w — akz]/{cketi^t) (as everywhere above, we use the braneh of (14) with the plus sign, therefore 
7^ = v^n2(l - /32)fc| + (n2 - (3^)k^/{nky/l - P^) > 0), and ai,2(k) are the amplitude operators of the TM^ and TE^ 
modes, respectively (the coefficients in front of (39)-(40) are to ensure that these operators have the dimension of 
\/m^). The unit vectors ei^2 are defined as ei = (k x zo)/|k x zo| and 62 = (k x ei)/|k x ei|. The Hamiltonian (38) 
can be now expressed as 



n = 



(27r)3 



^uj{k) [at(k)a,(k)+a,(k)at(k)], 



(41) 



where the index q — 1,2 labels the two main polarizations. 

The operators 01^2 (k) that we have introduced above for the two main polarizations must satisfy certain commuta- 
tion relations that should in the end lead to the commutation relation (31) for the quantum fields D(x, t) and B(x, t). 
We may write 



B(x,t) 



chut f d'^k^o^/k 



ei 



ai(k)e*('^-"-'^*) +al(k)t 



-«(k-x— cui) 





r d^kVk 


^2cst J 


' (27r)37o''' 



(42) 



D(x,t) = 



d^k 



2cAtt J (27r)37oV^ 
rf=^k7o 



(27^)3^/fc 



a2(k)e*(''"'-"*) + 4(k)e-'(''-^-"*) 
(ei X k) L(k)e*(''-"-"*) +al(k)e-'('""-'^*) 
(62 X k) fa2(k)e^('"^-'^*) +4(k)e-*(''-^-'^' 



+ 



(43) 



It can be verified that the commutation relation (31) follows from these formulas if the operators ai,2(k) satisfy the 
canonical commutation relations for annihilation and creation operators for bosons: 



[a,(k),a,(k')] = [at(k),at(k')] = 0, 
[a,(k),at(k')] = -[at(k),a,(k')] = {2n)%jSik-k'), 

where 5ij is Kronecker's delta. Indeed, with the help of the above formulas one may write 
t)B-(^'t)]--^ /^^k-VFf 7^(kxe,).(e,), 

[A(X,i),i^,(X ,i)J - 2 J (2^)3 J (2^)3 ^{ 

[ai{k),a\{k')y'-^-''-^'-'''-'^*+'^''^ + [al(k),ai(k')]e'(-'"^+*''"''+'^*-'"'*') + 



(44) 
(45) 



+ 



7o(k X 62)^(62)^ 
7^, 



([a2(k), 4(k')]e*('^-"-'''-"'-'^*+'^'*) + [4(k), a2(k')]e^(-'"-"+'"'-"'+"*-'^'*)) |, (46) 



where 7q = 7o(k') and ui' = ui{k'). Substituting (45) into (46) and taking the integral over k' one obtains 

d^k 



[A(x,i),B,(x',i)] = A^[k X (6161+6262)]^,^. (e'-^ ^— ') - e-^>^-(— ')) = 

_|_g.k.(x-x ) ^ _if^..^d^s{^ _ x'). (47) 



d 

dXn 



In this derivation we used the fact that the vectors 61,2 and k form a triplet of mutually orthogonal vectors, and, 
thus, k X (eiei + 6262) = k x 7, where I is the unity dyadic: (/)^^. = Sij. 

In a similar and simpler manner one can also check that the relations (44)-(45) ensure that all components of 
B(x,t), as well as all components of D(x, t), commute among themselves. 

Therefore, following [32, 33] we may conclude that the quantization of the electromagnetic field in a moving lossless 
nondispersive medium leads to the canonical result, with all the field operators and the Hamiltonian expressed in 
terms of the standard annihilation and creation operators of a bosonic field. 
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V. THE EXPRESSION FOR THE ZERO-POINT ENERGY 



The canonical diagonalized form of the Hamiltonian (41) allows for introduction of the particle number operator 
A/k,<j. By definition, the action of the number operator on a state results in the number of photons in this state: 
'^k.g = (^'k,g|-A/'k,g|*k,g)- Howcver, in Order for this to hold the states must be properly defined and normalized, so 
that (^'k,q|^'k,q) = 1- One way to achieve this is to discretize the k-vector space into cells of infinitesimal volumes Vk 

centered around the points k and require that in each state each cell contains an integral number of photons. 

Then, the number operator can be introduced as A/k.^ = Jy^-^^a^g(k')aq{k') = J-^^ ■j^^aj(k')aq(k) = 

/vk {2^)3 ^q{^)'^qO^')- The last two equahties are equivalent and hold because Vk is infinitesimal. Then, from the 
commutation relation (45) it follows that the Hamiltonian (41) may be expressed in terms of the number operator as 



(48) 



where the first sum is taken over all the cells in the wave vector space. In (48) we have labeled the frequency with an 
index q just to remind that the two main polarizations could in principle have different dispersion (which is the case 
of Kong's paper [33] where the medium at rest is assumed uniaxial). The term 

^ = EE^ (49) 

corresponds to the so-called zero-point energy, i.e., to the energy of the ground state of a quantum field. The sum 
(49) is wildly divergent and must be treated with a suitable renormalization procedure. It is known, however, that 
besides being divergent the zero-point energy £ may in some situations lead to physically observable phenomena, for 
instance, it plays a key role in the physics of the Casimir-Lifshitz forces. 

In a typical scenario in which one may observe a force due to the zero-point fluctuations of a quantum field, there 
exists a geometrical parameter, rf, that affects the modal dispersion and the density of quantum states of a system. 
Hence, this parameter, by virtue of (49), also affects the zero-point energy: £ = £{d). Any slow rate, quasi-stationary 
variations in this parameter result in variations in the amount of energy associated with the quantum fluctuations, 
which means that there appears a macroscopic force proportional to d£{d)/dd. 

For example, let us consider a layer of moving medium sandwiched in between two perfectly electrically conducting 
(PEC) mirrors positioned at a; = ±d/2. As before, we assume that the medium moves along the z-axis, so that the 
introdiiced mirrors do not interfere with the movement. It is evident that in this problem the modal spectrum is 
discrete in (to see this one has to complement the field equations (6)-(7) with the boundary conditions at the 
mirrors), while ky and kz form a continuous spectrum. Therefore, (49) may be written as 

£ _s-^sr' f '^^v ^(Wrf.fcy.fc.).g /rn^ 
- 2^ 2^ J -(2^2 2 ' 

q n ^ ' 

where £/L'^ has the meaning of the energy in the considered cavity per unit area of the mirrors; the infinite summation 
over n e Z skips n = for the TM^ modes. 

The frequencies that appear in the summation (50) may be understood as the eigenfrcquencies of a resonator 
formed by a layer of a moving medium and the mirrors. In general, for a pair of polarization sensitive, i.e., anisotropic 
mirrors (we will need this for the next section) the modes of such a resonator can be found by introducing 2x2 
reflection matrices (or, in other terms, planar dyadics) Ri^2{^- ky kz) of the mirrors and the complex propagation 
factor 7(w, ky, kz) = —ikx{oj, ky, kz) of the waves that travel in between the mirrors. Then, in terms of these quantities 
the characteristic equation for the modes in the cavity is readily obtained as 



V{uj,ky,kz,d) = Aet I^^) - Rx{w,ky,kz) ■ R2{oJ,ky,kz)e-'^'<^'^^^y'''^^'^ =0, (51) 

where 7(^) is the planar unity dyadic. The characteristic equation (51) for the case of the ideally conducting mirrors 

reduces to (1 — e"^'''')^ = with the obvious solution = nn/d that appears in (50). 

In what follows we are going to use the principle of argument to replace the discrete summation over the resonant 
frequencies in (50) by an integration in the complex plane of oj. Indeed, if a function /(w) is analytic and has 
roots in a closed region G with the boundary dG, then the sum over its roots in this region can be found as 
^u)k = (27ri)^^ §QQU!d\ogf{Lu). There is, however, a subtle difficulty when applying this principle to the function of 
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the characteristic equation (51), because ky, kz, d) may have poles and branch points. The poles may appear at 

the points where the reflection coefficients i?i,2(w, ky, kz) have resonances: \Ri,2\ — ^ cjo, and, thus, they correspond to 
surface waves that may exist at the boundaries between two different media. Respectively, the branch points appear 
at the frequencies where ^{uj, ky, kz) = 0, i.e., at the points where the propagating waves transition into the evanescent 
ones. The main difficulty is with the branch points, as the poles of the reflection coefficients do not depend on d and 
only add a constant to the sum representing the zero-point energy, i.e., they only shift the origin of the zero point 
energy which is irrelevant. 

However, it is possible to rewrite the characteristic equation in a form that does not have branch points and is 
mcromorphic in the complex plane of w (see Appendix A) . When such a form of the characteristic equation is used 
(here we use the same symbol V to denote the function of this characteristic equation), the summation over the 
discrete frequencies for a given pair of ky, kz in (50) can be formally replaced with an integration over a path C in 
the complex plane of lo that encircles the roots of the characteristic equation: 

= ^ / .d\ogV{u:,ky,kz,d). (52) 

When /9^n^ < 1 the roots we are interested in lie on the positive half of the real axis (see Section II), therefore we 
can choose the path C so that it follows the imaginary axis from +ioo to —ioo and then closes in the right half of the 
w-plane with a semicircle Coo of an infinite radius. 

It should be well understood at this point that the integral (52) diverges, as well as the original series (50) does. 
Nevertheless, one may find a way to regularize (52) by dropping distant-independent infinite terms in the integration 
(52), as explained in Appendix A. Doing this, the regular part of the zero-point energy, or, in other terms, the Casimir 
interaction energy at zero temperature, 5£, can be expressed with an integral over the imaginary axis only: 

+ioo 

5£ h f dkydkz f , , ,s ^ f dkydkz 

T^=-^^Jl2^ J -d^os^i-M,d) = -J 

— ioo 

+ 00 +00 

X j \ogV{ii,ky,kz,d)di=^j j \ogV{ii,ky,kz,d)di, (53) 

-oo 

where we replaced the integration variable by w = and integrated by parts once. The last equality holds due to 
the symmetry with respect to a simultaneous change of signs of ^ and k^'- T>{—i£^, ky, —kz) = T>{i^, ky, kz). One may 
recognize in (53) the so-called generalized Lifshitz formula for the Casimir energy at zero temperature. 

As we are not using a covariant formulation of electrodynamics in this paper, the relativistic covariance of the 
obtained result (53) requires an additional discussion. Some implications of special relativity on the reflection matrices 

i?i,2 are outlined in Appendix B. In particular, it can be verified that if there exists a reference frame in which the 
moving media are at rest, then (53) reduces to the classic Dzyaloshinski-Lifshitz result [13] for the Casimir energy of 
stationary magnetodielectric slabs. This is, of course, just a consequence of the material relation transformations (3)- 
(5). 

Above the threshold of the Cherenkov radiation, i.e., when P^'n? > 1, the dispersion relation (14) may result in 

negative frequencies irrespectively of which branch of (14) is selected. By virtue of (48) this leads to the appearance 
of negative quanta in the range of wave vectors that belong to the Cherenkov cone. These quanta are responsible for 
a potential instability in a medium that moves with a velocity higher than the phase velocity in the same medium at 
rest. Indeed, any stationary (in the laboratory frame) object that perturbs the electromagnetic field will radiate in 
such quickly moving medium. Because of this unavoidable instability, in the following sections we restrict our analysis 
only by the case when ^^n^ < 1. 

VI. CASIMIR ENERGY AND FORCE IN LAYERED MOVING MEDIA 

In this section we will consider the Casimir energy and force that result from the zero-point fluctuations in a layered 

moving medium, namely, in a configuration analogous to the canonical problem solved by Lifshitz [2, 13]. Thus, we 
consider a structure composed of a moving layer of finite thickness d sandwiched in between two other semiinfinite 
moving layers (Fig. 1). The velocities are assumed uniform within the layers and collinear with the z-axis which is 
parallel to the layer interfaces. Such a structure may be understood as a simplified model of a nonuniformly moving 
fluid, in which the width of the transition regions where the velocity changes continuously is assumed small compared 
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Ri 



(ii) 
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FIG. 1: (Color online) A layer of moving magnetodielectric (iii) of thickness d sandwiched in between two moving semiinfinite 
magnetodielectric layers (i) and (ii). The three layers slide along the same line which is parallel to the interfaces of the layers. 
The magnitudes and the signs of the velocities i'i,2,3 are arbitrary. The Casimir-Lifshitz force is calculated from the reflection 

matrices 7?i,2 defined at the interfaces of the layers. 



to the thickness of the layers. In other words, we neglect all friction effects that may exist at the boundaries of the 
moving layers. In practice such situation may be achievable, for example, in certain phases of liquid helium at very 
low temperatures or within metamaterial layers where the velocity is merely a structural parameter (i.e., when there 
is no real movement). 

It is evident that one may always choose a reference frame in which the middle layer is at rest. However, we prefer 
not to impose such a restriction, ensuring in this way a straightforward generalization of our results to the case of 
multiple moving layers. As follows from the treatment of Section V, in order to obtain the Casimir energy of this 
system, one must first solve for the reflection matrix at an interface of two moving layers. Although there are some 
results available in the literature (see, e.g., [40] and references therein), they are typically given in a form unsuitable 
for our purposes (e.g., one of the layers is assumed to be vacuum) therefore, in Appendix B we derive the necessary 
expressions for the components of the reflection matrices 

Ri.2 = ( ^hh j (54) 

that are defined in terms of the z-componcnts of the electric and magnetic fields. With these expressions at hand, the 
Casimir interaction energy in the canonical triple-layer structure is given by (53) where the matrices i?i.2 correspond 
to the two interfaces of the middle layer. From the expressions derived in Appendix B, it is also seen that the reflection 
matrices are invariant under a simultaneous change of signs of lj and ky/. we used this property when obtaining the 
expression for the zero point energy (53). 

Next, the Casimir force component normal to the interfaces is found by differentiating (53) with respect to the 
thickness of the middle layer d (we use the convention that a positive force corresponds to attraction): 

oo 


h ^ fdky dk, f 2A„7e-2T"^ 

n-L g 

where Ai^2 are the eigenvalues of the matrix Ri • i?2- It can be shown that the same expression for the Casimir force 
must also hold in the case of dispersive material parameters which is discussed in Appendix C. 

Because the dispersion equation for the waves in (lossless and nondispersive) moving media is not symmetric with 
respect to the change of sign of the frequency u, it is not anymore a function of o;^ as in (lossless and nondisper- 
sive) reciprocal media. Due to this asymmetry the reflection matrices i?i.2 are in general complex at the imaginary 
frequencies uj = while the respective matrices in reciprocal media are always real under the same circumstances. 

Therefore, in general, the eigenvalues Ai,2 of the matrix i?i • R2 are also complex. In Section VII we will, however, 
show that the expression for the Casimir force always results in real numbers, due to the symmetry of the integrand 
of (55). 
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FIG. 2: (Color online) The absolute values of the reflection coefficients -Rabs — l-R'^'^l and -Rljbs — (the respective plots 

coincide and are shown with a single surface) at the real frequencies as functions of the normalized wavenumbers kz/ko and 
ky/ko at an interface of a stationary medium with Sr ~ 2, fii = 1 and the same medium moving with velocity v = 0.6c. The 
plotted surface is colored proportionally to the reflection amplitude, as indicated in the color bar on the left. 



To simplify the integral (55) further we introduce new dimensionless variables Hy = cky/^, Kz = ck^/^, v — cj^j^., 
and C = ^c^/c, in which (55) becomes 

_ he f dKy dKz f 2t/A„C'e-^-g 

L2 - 2^d4 J i2nY J 1-A„e-2K ^^^^ 

n-l (, 

One may notice that both A„ and v do not depend on Q (they depend only on the relative wavenumbers Ky and 
because the material parameters are assumed nondispersive), therefore, by substituting Q = t/{2i') we obtain 



OO 



Fc he /■ dny dKz Xn f t'^e * dt 3hc f dKy dnz Li^^Xn) , , 

1? " 16^ J (27r)2 ^ J A„e-* " 8^ J "(2^ ^ ' 

where the integral over t results in the polylogarithm Li4(z) = X^^i /n'^. 

Thus, the Casimir force in layers of moving (nondispersive) media has the same dependence on the distance as the 
Casimir force between two ideally conducting plates in vacuum. It is also seen that the value and the sign of the force 

(57) are determined by the the eigenvalues Ai^2 of the matrix Ri ■ i?2, which in turn depend on the relative velocities 
of the layers. In the next section we will study numerically this dependence and will demonstrate that the Casimir 
forces in moving media may be repulsive under certain conditions. 



VII. NUMERICAL EXAMPLES AND DISCUSSION 



In this section the expression for the Casimir force (57) is analyzed numerically. It is convenient to start from 
discussing some properties of the reflection coefficients (54). First of all, we would like to remind that the elements 
of the reflection matrix (54) are defined in terms of just a single component of the electric and magnetic field vectors 
(see Appendix B). Therefore, in general, their values differ significantly from the classic reflection coefficients into 
CO- and cross-polarized TM and TE waves (the cases when (54) reduces to the classic formulas are mentioned in 
Appendix B). For example, the magnitudes of the cross-components R^^ and R^'^ in our definition may exceed unity 
when the characteristic impedance of the layers is different from the free-space impedance 770. 

At real frequencies the elements of the refiection matrix (54) behave as shown in Figs. 2-3. In these figures we plot the 
absolute values of the reflection coefficients at an interface of a stationary medium with the relative parameters = 2, 
/ij. = 1 and the same medium moving with velocity v = 0.6c along the z-axis as functions of the relative wavenumbers 
ky/ko Bind kz/ko (where ko = uj/c) of an incident wave (the wave is incident from the side of the stationary layer). 
The cross-components of the reflection matrix plotted in Fig. 3 are normalized as indicated in the flgure caption. In 
these figures only the propagating waves are considered, i.e., the waves with {ky/ko)'^ + (kz/ko)^ < Er^r- 

As one may notice, the elements of the reflection matrix demonstrate a strongly nonreciprocal behavior: the 
reflection is different for the incident waves with positive and negative kz- It is also noticeable that the reflection is 
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FIG. 4: (Color online) Real and imaginary parts of the eigenvalues Ai,2 of the matrix Ri ■ R2 at imaginary frequencies as 
functions of the normalized wavenumbers and Hy (the plots for Ai and A2 coincide and are shown with a single surface). 
There are three layers of the same medium with £r = 2, /Xr = 1. The middle layer is stationary and the two outer layers move 
with the velocity v = 0.6c along the positive direction of the z-axis. 



rather low overall because the parameters of the layers are chosen so that there would be no reflection if there were no 
movement. However, the grazing waves reflect strongly, as well as the waves with kz/ko > {y/StUt — a)/ y/solM) ~ 1.09. 
The latter is due to the fact that the waves with k^ greater than the mentioned limit are evanescent in the moving 
layer, as can be easily seen from the dispersion equation. 

The behavior at the imaginary frequencies is better understood from the eigenvalues Ai^2 of the matrix Ri ■ R2 
written for the complete structure composed of the three moving layers. Accordingly to (57), these eigenvalues 
determine the magnitude and the sign of the Casimir force. The plots of the eigenvalues are given in Fig. 4 for the 
case when the outer layers move in the same direction with velocity v = 0.6c, and in Fig. 5 for the case when the two 
outer layers move with the same speed, but in the opposite directions. The middle layer is stationary in both cases. 

In the case when the two outer layers move in the same direction with the same velocity (Fig. 4) the two eigenvalues 
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FIG. 5: (Color online) The eigenvalues Ai,2 at imaginary frequencies as functions of the normalized wavenumbers Kz and Hy 
for the case when the outer layers move in opposite directions (the eigenvalues are purely real in this scenario). The absolute 
value of the velocity and the other parameters are the same as in Fig. 4. 





FIG. 6: (Color online) The real and the imaginary parts of the integrand of Eq. (57) as functions of the normalized wavenumbers 
Kz and Ky in the same scenario as in Fig. 4. 



of the matrix Ri ■ i?2 coincide. The eigenvalues are complex in this case, with the real part concentrated mostly in 
the negative half space, and the imaginary part changing sign when changes sign, which is a consequence of the 

fact that (^Ri^2ii^, ky, kz)j = Ri^2i—iS,,ky,kz) ~ Ri^2{iS.,ky,—kz) when ^, ky and kz are real. 

When substituted into the integral (57) the dominating negative real parts of the eigenvalues result in a negative 
Casimir force, which corresponds to a repulsion. The contribution of the imaginary part vanishes due the symmetry 
of the integrand. To further illustrate this, in Fig. 6 we plot the integrand of (57) as a function of the normalized 
wavenumbers Kz and Hy. As is seen, only a small area of the {Kz,Ky) plane contributes to the integral, with the 
negative values of the integrand on the periphery of this area clearly outweighing the positive values seen at the 
middle. 

When the two outer layers move in the opposite directions with the same absolute speed (Fig. 5) the eigenvalues 

of the matrix i?i • i?2 are both real and positive (in a less symmetric scenario when the absolute velocities of the 
two layers differ there also appears a non-zero imaginary part). Thus, this case results in attraction between the two 
moving layers, as clearly seen from the plot of the integrand of (57) in Fig. 7. This agrees with findings of Ref. [3], 
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FIG. 7: (Color online) The integrand of Eq. (57) as a function of the normalized wavenumbers and Ky in the same scenario 
as in Fig. 5 (the integrand is purely real in this scenario). 
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FIG. 8: (Color online) The magnitude of the attractive and repulsive Casimir-Lifshitz forces in the triple-layered structures 
with = 1 and Sr indicated in the plot, as functions of the relative velocity v/c of the outer layers (logarithmic scale). The 
force is normalized to the Casimir force between two perfect electric conductors (PEC) separated by the same distance as the 
thickness of the middle layer. The arrows in the plot indicate the directions of the movement of the outer layers that result in 
attraction and in repulsion. 

where only this type of relative movement of dielectric slabs (separated by vacuum) was considered. 

To further study the attraction and repulsion phenomena in moving layers we have calculated the velocity depen- 
dence of the attractive and repulsive Casimir-Lifshitz forces in the two scenario considered above. The results are 
represented in Fig. 8. In this figure we plot the magnitude of the force \Fc\ normalized to the attractive Casimir force 
in a system of two ideally conducting plates i^pec = 7r^/ic/(240(i'*), where d equals the thickness of the middle (sta- 
tionary) layer (as we noticed in Section VI the dependence of the force on distance in layers of moving nondispersive 
dielectrics is the same as in Casimir's canonical structure). Fig. 8 also demonstrates the dependence of the force on 
the value of the dielectric constant. 

One can see that at low velocities the repulsive and attractive forces in the two scenaria of the relative movement 
of the outer layers are close to each other, while at larger speeds the attraction is stronger than the repulsion. The 
double logarithmic scale of Fig. 8 indicates that at small velocities both forces are proportional to (u/c)^, thus, the 
effect reported in this paper has the same order as most of the relativistic effects. Quite naturally, the effects are 
more pronounced in media with higher permittivity. 

In the last numerical example we calculate the attractive force between a stationary and a moving dielectric 
separated by a vacuum and compare it with the same force derived in Ref. [3] with an independent Green tensor- 
based approach. The results of this comparison can be seen in Fig. 9, where we plot AF — F{v) — F{y — 0) which is 
an addition to the force that appears because of the relative movement of the layers. 

From these calculations we conclude that up to the accuracy of numerical integration (which is a triple integration 
in the case of Ref. [3]) the results of the two independent approaches expressed by Eq. (57) of the present paper 
and Eq. (40) of Ref. [3] are in excellent agreement. The same reference contains also an expression for the leading 
0(/3^) term of the velocity-dependent correction to the Lifshitz force (Eq. (42) of Ref. [3]). However, one must be 
accurate when making a comparison against this result, because the (velocity-dependent) addends A'^\, and 
seem to appear there not expanded in powers of /3. A plot of the explicit /3^-proportional term [42] of the mentioned 
expression is shown in Fig. 9 with a blue dashed line which docs not match the exact result at low velocities. 
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FIG. 9: (Color online) The additional attractive force AF = F(v) — F{v = 0) exerted on a dielectric with the relative 
permittivity Ei moving with the relative velocity v/c nearby a stationary dielectric of the same permittivity, for the three 
different values of the relative permittivity: 2, 4, and 8. The dielectrics are separated by a vacuum gap. The force is normalized 
to the Casimir force between two stationary PEC plates separated by the same gap. The brown solid lines: the force calculated 
with the theory of the present paper [Eq. (57)]. The blue dots: the same force calculated from Eq. (40) of Ref. [3]. The blue 
dashed line: the plot of the /3^-proportional term of Eq. (42) of Ref. [3] for Si — 2. 

Although it is out of the scope of this paper, the observed agreement suggests that calculations of Ref. [3] are 
applicable to the geometries that can be considered as effectively closed ones (which is also the case of this paper) 
in which the pertinent difficulty with the branch cuts pointed out in Refs. [4, 6] can be treated in a manner similar 
to what we have done in Appendix A. Indeed, in this work we have shown that the branch points of the reflection 
coefficients of moving layers are irrelevant in such geometries. 

VIII. CONCLUSIONS 

In this paper we have considered the forces due to quantum-mechanical fluctuations of the electromagnetic field 
in layered moving media. We have demonstrated that rapid relative movements of neighboring layers in a dielectric 
(e.g., in a nonuniform fluid flow) may result in both attractive and repulsive interactions between the layers. 

Although in the present study we have made an emphasis on the Casimir-Lifshitz forces resulting from relativistic 
movement of material layers, the results of this paper apply also (at least, qualitatively) to a class of bianisotropic 
metamaterials called moving media. Thus, we may conclude that a speciflc type of nonreciprocal magnetoelectric 
interaction in bianisotropic composites may also result in repulsive Casimir-Lifshitz interactions. There have been 
previous attempts to realize Casimir repulsion in metamaterials with the help of reciprocal magnetoelectric interaction 
(e.g., chirality). However, it was recently shown [25-27] that the causality and passivity preclude Casimir repulsion 
in reciprocal metamaterials. 

The Casimir-Lifshitz interactions studied in this paper may be of importance in areas of physics involving rapid 
movements of matter, as well as in the phenomcnological quantum electrodynamics of nonreciprocal materials. 
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Appendix A: 

The problem of branch points in the context of Casimir's energy calculation dates back to 70's of the last century. 
Some of the main ideas of the approach that we are going to use in this Appendix have been borrowed from Ref. [41]. 

Instead of considering an initially open structure, we start with the situation in which the moving layers are bounded 
by PEC walls, as depicted in Fig. 10(a). As is seen, there are two PEC-backed layers of media (i) and (ii) that can 
move in a background filled with medium (iii) that is in turn terminated by two PEC walls at a; = and x = L. We 
assume that d ^ bi 2, and &i_2 ^ L. When the PEC-backed layers (i) and (ii) move, their thicknesses bi 2, as well as 
the total size of the structure L — ci + bi + d + b2 + C2, remain fixed. 
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FIG. 10: (Color online) (a) The PEC-backed structure used in calculation of the distant-dependent part of the zero-point 
energy, (b) The integration path C in the complex plane of uj. 



In this structure, the three regions 0<x<ci, ci < x < L ~ C2, and L ~ C2 < x < L are electromagnetically 
screened from each other. Therefore, the characteristic equation for the whole structure is a product of the equations 
for the three regions: 



(Al) 



In the middle of (Al) one can recognize the term that has the form (51); we have also made explicit the dependence 
of terms of (Al) on the propagation factor in the background medium 7. 

The reflection coefficients i?i,2 have the following important property: 



(A2) 



which can be seen from the fact that the reflection dyadics can be expressed as i?i,2(7) = — 

1 — 



Z 



1,2 



I - z 



1,2 • 



I + -^1,2 • Ywil) , where Zi^2 are the dyadic input impedances of the 

PEC-backed layers which are meromorphic in the whole complex plane of lo and independent of 7, and Kw(7) is the 

dyadic wave admittance of the middle layer that is such that Y^{—^) = —Y^{'^). It should be noted here that while 
the reflection matrix of an open half space has the same property (A2), the input impedance of such a space is not a 

meromorphic function of uj (see Rcf. [41]). Thus, we may conclude that the branch points of Ri^2 coincide with the 
branch points of that are at the frequencies where 7(w, ky, k^) = 0. 
Using the above property we may express I?(— 7) in terms of I'(7): 



I?(-7) = (l-e 



+27C1 



det <^ / - i?i 



h27d 



9 f)('y\p'^l<.I--bl-b.2) 



det 



{i?i-i?2} 



(A3) 



where Ri^2 = -^1,2(7)- From (A3) one can see that the roots of function 2?(— 7) in w include, in general, all the roots 
of 'D{'~f). Thus, we may construct a function 



J-(7) = ^(7)^(-7) 



det|i?i •i?2} 



det \Ri- R2 



} 



(A4) 
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where V(-i) = det |/- iii(7) • R2{l)(i'~'^^'^'^ has the form (51). The function F{-f) has all the roots of V{-f) (with 

a difference that simple roots of V become roots of second order in and is even in 7. The latter makes ^-"(7) a 
meromorphic function of oj. 

Therefore, we may apply the principle of argument (as explained in Section V) to this function instead of applying 
it directly to X'(7). The integral over the respective path (see Fig. 10(b)) in the complex plane of w reads in this case 



/wdlogJ"= / wd\ogV+^. I wdlog[(l-e-2^'=0''(l 
Am J J Am J ^ ^ 



Cab Cab 



I r pil(L-bi-b2) 

+ 73^ / ^d\og ' ^ (A5) 

det|iii-i?2| 



Cab 



where Cab is an open path that is obtained from C by introducing a cut at the point where the semicircle crosses 
the real axis. Such a cut is necessary because the expressions under the integrals on the right hand side of (A5) are 
not meromorphic in uj. 

Physically, the integral (A5) represents a part of the zero point energy that is due to the modal frequencies which 

are the roots of (Al) that have been encircled by the path C. Because we are interested only in the variation of the 
zero point energy with the separation d between the two moving slabs, we may drop the last addend on the right 
hand side of (A5) as it is independent of d. The second addend can be made arbitrary small when ci,2 — >■ 00 due to 
the nonvanishing positive real part of 7. Thus, the distant-dependent part of the integral (A5) is given by 

+iR 

-dlogV = -l-. j .dlogV+±. J .dlogV, (A6) 

Cab —iR Caa'^C^/^ 

where the first integral on the right hand side is taken over a path that lies on the imaginary axis and the second 
integral is over the two halves of the semicircle. 

The integral (A6) depends on the thicknesses of the slabs bi 2 and the slab separation d. Now we let 61 2 — >■ 00 
in (A6) (when taking this limit, we assume that still L ^ ^1,2)- In this limit, due to the nonvanishing imaginary part 

of oj under the integrals on the right hand side of (A6), the reflection coefficients of the PEC-backed layers i?i^2 will 
tend to the respective reflection coefficients of open half spaces (which are derived in Appendix B). 

The last step of the derivation is to let the radius of the semicircle tend to infinity: i? — >■ 00. In this limit, which 
corresponds to infinitely high frequencies, all dispersive materials (including the materials with very weak dispersion 

that we consider in this paper) become transparent. Therefore, i?i_2 — > under the integral over the semicircle, and 
this integral vanishes. This leads to the expression (53) for the interaction part of the zero-point energy. 



Appendix B: 

Let us consider an interface in a pair of layers. Without any loss of generality we let the interface be at a; with 
the a;-axis orthogonal to the interface. At the interface the tangential components of the electric and magnetic fields 
of the two main polarizations are given by Eqs. (lO)-(ll): 

^ w^EtAtt - (kz + wa)2 
^ kyih+^a) (TM,) (Bl) 

^ w^etMt - {kz + oja)"^ 



ij'^etfJ-t - (kz + ujaY 



ky {kz 



w^£tMt - {kz + wa)' 



;Hz 



Ez=0, 



(TE.) 



(B2) 



19 



where we have replaced k"^ in the denominator with an equivalent expression that follows from Eqs. (8)^(9). These 
relations hold at both sides of the interface, with the material parameters et, fJ-t, and a, and the wave vector components 
taken at the respective sides. 

In the following we are going to formulate and solve a plane wave reflection problem at an interface of two moving 
media. To simplify writing we introduce the following notations 

^ ^ kyjkz + kpa/y/eofio) 

^ feofex(gt/£o) 

kletHt/{eoiip) - (k^ + fcoa/\/eoMo)^ ' 



oH kok^{pt/fJ-o) 



koStl^t/ (eoMo) - {kz + koa/ y/e^f ' 



(B5) 



where fco = w^eoMo- Then, with these notations at hand we consider a TM^ wave of unit amplitude incident from 
the region a; < and write the fields in this region (the factor 6*^*'''^+*^'^) common at both sides of the interface is 
dropped) as 

Ez = e**^-''^ + Ae-''=»''^, (B6) 
Hz = Be-'''-^'', (B7) 
Hy = -,,-i/3f e^'^^"" + ?7o"'/3f ^e-''^^''^ + alBe-*'=^''^ (B8) 
Ey = aie*'=^"=" + aiAe-''^^"^ - r?o/3f 5e-**=^"^, (B9) 



and in the region a; > as 



E, = (7e'^»"^, (BIO) 

ff^ = De'''-^'^, (Bll) 

Hy = -rt^^p^Ce^-'"'' + aiDe''"-^", (B12) 

Ey = aaCe^'^^''^ + r/o/Sf £>e''=«''^ (B13) 



where A, B, C, and D are y et unknown wave amplitudes of the two reflected and the two transmitted waves, 
respectively, and rjo = \J jJ^a/sQ. As one can see, we take into account the fact that a TM;^ incident wave may produce 
in general both polarizations in the reflected and transmitted fields. 

Equating the tangential components of the electric and magnetic fields at both sides of the intcirfacc; at x one 
obtains a system of four equations for the four unknown wave amplitudes. Solving this system for A and B (i.e., for 
the reflected waves) we find 

(ai-a2)^-(/3f-/3|^)(/3f + 

(«i - a^f + (/3f + /3f )(/3f + /3f ) ' ^ ' 

(a,-a2)2 + (/?f + /3f)(/3f +/3f)- ^ 

The case of a TE^ incident wave can be considered in a completely analogous manner. Below we give just the final 
result for the amplitudes of the refiected waves: 

w_ (ai-a2)^-(/?f -/?f)(/3f + /3f) 

(ai-a2)2 + (^f + /3f)(/3f +/3f)' ^ 
o/ 2?7o(qi - Q;2)/3f 

(ar-a2)2 + (/3f + /3f)(/3f + ^ 

Thus, we may introduce the following refiection matrix written in terms of the ^-components of the fields: 

Ef \^( A rj^'B'\ ( Er \ = fR'" R'""] ( Er \ (^..^s 
rjoHf UoB A' ; • [rioH'r " U'^° i?"V ' VloH'r ' ^ ' 
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As can be verified, the elements of the reflection matrix reduce to the standard Fresnel reflection coefficients of 
the P- and S-polarized waves in the special case oi a = 0, ky = 0, for which R'^^ = R^'^ = 0, = Rp = 

(sik^x ^ - £2^1^^) / (eik^^ + £2fci^^) and R^^ = -Rg = [jiik^x '' - At2fci^^^ / {jJ'ik'i^ + iJ,2ki^^^ , and also in the case of 

a = 0, kz = 0, for which R^^ = R^^ = 0, R"^ = Rg, and R^^ = —Rp. In the general case, the standard reflection 
matrix defined in terms of the tangential components of the electric field can be obtained from the matrix (B18) with 
the following similarity transformation: 



(B19) 



As mentioned in Section VI, the reflection coefficients (B18) are in general complex, even at purely imaginary 
frequencies. The complexity of the reflection matrix (B18) at imaginary frequencies is an unusual property that by 
itself deserves a separate study. Here we will only briefly outline the main reason behind this complexity. Indeed, 
from a physical point of view, the reflection at imaginary frequencies uj — iS^ can be understood as the response to 
an incident wave that has the time dependence of the form e^*, i.e., to a signal that grows exponentially with time. 
Let us now consider an interface between a vacuum at a; < and a moving medium at x > 0, and let us assume that 
there is a plane wave with such time dependence impinging on the interface from the side of the the vacuum. We 
set up the same coordinate system as above so that the movement is along the ^-axis. In this coordinate system the 
incident wave of, for instance, the TM^ polarization can be written as 

E'^" = i;oe*('=''^+'='^)e«*-^^, (B20) 

where ky and kz are the real propagation factors in the interface plane, and 7 = —ikx = ^^^eoMo + ky + k^ > is 
the solution of the vacuum dispersion equation at imaginary frequencies. As we arc interested only in an illustration, 
we let fcy = in (B20), so that the TM^ wave becomes the standard TM wave with respect to the plane of incidence. 

A vacuum is invariant under the Lorentz transformations (see Section II), as are the components of the electromag- 
netic fields parallel to the velocity vector (the z-componcnts) , therefore to solve for the reflection coefficient we may 
switch to the comoving frame in which the reflection coefficient is simply 

where 

^, ^ i + ikzv ^, ^ K - iiv/c^ 

are the imaginary frequency and the z-component of the wave vector transformed to the comoving frame. Substitut- 
ing (B22) into (B21) we obtain after some manipulation 



V(n2-l)(e/c)2 +72+6,7' 

where 7 = ^/^'^Eqij.o + k"^ and = e/eo. As is readily seen, i?®'' is in general complex when ^ 1 and v 0, and 
this complexity is due to the fact that the relative movement of the layers intermixes the imaginary frequencies with 
the real-valued wavenumbers by the virtue of the Lorentz transformations. It is easy to check that the result (B23) 
is a particular case of more general formulas (B14)-(B18). 

Conversely, one may verify that if there exists a reference frame at which the moving matter is at rest, then under a 
transformation of the form (B22) the complex propagation factor ky, k^) and the reflection matrix (B18) reduce 
to the respective expressions in stationary magnetodielectrics. Additionally, when such a transformation is applied 
to the integrand of (53), one may notice that the integration element dkydk^d^ is preserved, because the Jacobian 
of the transformation (B22) equals unity: d{^,kz)/d{^' ,k'^) = 1. Therefore, the Casimir force per unity of area (the 
Casimir pressure) given by (53) is the same in all reference frames that move parallel to the layers, provided that the 
velocities of the layers arc transformed accordingly to the rclativistic velocity addition law. Such an invariancc of the 
Casimir pressure (53) is not surprising, as physically the pressure exerted on the moving layers is related with the 
component of the photon momenta that is perpendicular to the direction of the movement, and this component is 
preserved under the Lorentz transformation. Thus, we may conclude that our formulation extends the known theory 
of Casimir-Lifshitz forces in dielectric layers in a way fully consistent with special relativity. 
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Appendix C: 



In this appendix we discuss how the results obtained for non-dispersive moving media may be generahzed to include 

the effects of frequency dispersion in the effective material parameters. 

Let us consider an isotropic dispersive magnetodielectric described by the following material relations in its proper 
frame: 



oo 

T)'{^,t')=so j £,(T')E'(x',t'-r')dr', 



(CI) 



CSO 

B'(x',t') = MO y Mr(r')H'(x',t' - r') dr' , (C2) 



where erC^') and Hi{t') are the dielectric and magnetic response functions. 

In the proper frame which is co-moving with the medium, the field components orthogonal to v can be expressed 
through the same components in the stationary laboratory frame as 

E;=7L(Et+vxBt), H;=7L(Ht-vxDt), (C3) 



D;=7L(Dt + ^vxHt), B[ = 7L(Bt - X EO, (C4) 

where v is the medium velocity (along Oz) and 7l = l/-\/l — v'^/c^. Substituting (C3)-(C4) into (C1)-(C2) one 
obtains 



oo 

Dt + X Ht = £r(r') [Et(^(^',t'),i(^',t' - r')) + v x B^{z{z' ,if),t{z' ,t' - r'))] dr', 



(C5) 



OO 

Bt-^vxEt=Moy tir{r')[llt{z{z',if),t{z',i/ -T'))-vx-DMz',t'),t{z',i^ -r'))]dT', (C6) 



where z = z{z', t') = ^\,{z' -\- vt'), t = t{z', t') = 7L(t' + vz'/cP). Prom here, 

OO 

Dt + X Ht = eo J £r(T') [E^{z - jj^vt', t - 7lt') + v x Bt(^ - 7LVr', t - 7lt')] dr' , (C7) 



Bt - X Et = Mo j Mr(T') [Ht(^ - 7Lt^T',i - 7lt') - V X Dt(2 - li^vr' ,t - 7lt')] dr' . (C8) 



In order to obtain the constitutive relations in the laboratory frame, one has to solve the system of integral equations 
(C7)-(C8) to express D and B in terms of E and H. It is evident that, in general, the above system may not result 
in a simple proportionality relation between the flux and field vectors. However, for plane waves this system is easily 
solvable and results in relations (l)-(5) of Section II with e = e{oj'), n = m(w'), and n = n(w'), where uj' = ji,{uj — kzv) 
is the angular frequency in the proper frame of the moving nicidium. As this frequency depends on the wavenumber in 
the laboratory frame, the relations (l)-(5) with the modified parameters readily describe a spatially nonlocal medium, 
as was mentioned in Introduction. 

One may also verify that such modification does not affect the frequency domain treatment of Section II. The 
equations (8)-(13) written for the plane waves in a moving nondispersive magnetodielectric hold also in the case of 
dispersive moving media if the parameters et = £t(w'), A*t = A'tC'*^'); ^ — are understood as w'-dependent. 
Eq. (14) becomes a transcendental equation in the dispersive case. The important symmetry of Eqs. (8)-(13) with 
respect to the simultaneous change of signs of ui and k discussed in Section II is preserved in the dispersive case. 
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because k) = — cj'(— — k) and e{—uj') = s*{uj'), w') = Thus, the generahzation of the classical part 

of this study to the dispersive case is trivial. 

The quantum-theoretical part of this paper is based on the expressions (27)-(28) and (38) for the Haniiltonian 
of the free electromagnetic field. As has been mentioned in Section III, these expressions are physically understood 
as summations over the energies of all possible modes in a modal expansion of the electromagnetic field. Therefore, 
it is only natural that the same expressions must also hold in the case of frequency dispersive material parameters, 
provided that the basic relations for the energy w and the momentum p of a photon in a dispersive medium remain 
the same as in a vacuum: w = Hco, p = Hk, w/p = uj/k = Uph- Hence, one must also expect the diagonalized form 
of the Hamiltonian (41) to be valid in the dispersive case, in which the modal frequencies a;(k) are found from the 
transcendental equation (14) that must take into account the material dispersion. 

Therefore, the expressions for the interaction part of the zero-point energy (53) and the Casimir force (55)-(56) 
must also hold in the dispersive case. 
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